Jaynes-Cummings model: Counter rotating effect on the vacuum Rabi splitting and 

atom-cavity dynamics 
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The effect of the counter-rotating terms in the Jaynes-Cummings model is investigated with an 
extended coherent-state approach. The counter-rotating terms greatly modify the vacuum Rabi 
splitting. Two peaks with different heights in the weak coupling regime and more than two peaks 
in the intermediate coupling regime are predicted. The collapses and revivals in the evolution of 
the atomic population inversion disappear in the intermediate coupling regime, but reappear in the 
strong coupling regime. This reappearance is similar to that under the rotating-wave approximation, 
attributed to the summation of periodic cosine functions of the evolution. 
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The matter-light interaction is a fundamental one in 
optical physics. The simplest model is the Jaynes- 
Cummings (JC) model 1], which describes the interac- 
tion between a two-level atom and a single mode of the 
quantized electromagnetic field. In the semiclassical de- 
scription, the evolution of the atomic population inver- 
sion shows a regular oscillation behavior. In the quantum 
mechanics description with the rotating-wave approxi- 
mation(RWA), JC model is exactly solveableQ. By the 
exact solution, the revival and collapse of the atomic pop- 
ulation inversion has been predicted by Eberly et alQ- 
5j, and later confirmed experimentally Q. The other im- 
portant quantum mechanical feature is the vacuum Rabi 
splitting 0, [H , which is equal to 2g (g the atom-cavity 
coupling strength) in the JC model with the RWA. 

The RWA is valid in the weak coupling regimes with 
small detuning in the study of the JC model, because 
the contribution of the counter-rotating terms (CRTs) 
to the evolution of the system is quite small. In the 
coupling is not weak, the CRTs, which allow the simul- 
taneous creation and annihilation of an excitation quan- 
tum in both atom and cavity mode, must be taken into 
account. The CRTs would modify the above two inter- 
esting quantum phenomena, i.e. the revivals and col- 
lapses of the atomic population inversion and the vac- 
uum Rabi splitting. A perturbation by path-integral 
method has been applied to study the contribution of 
counter-rotating terms in quantum Rabi oscillation in 
intermediate-coupling regime [{j. The effect of the CRTs 
in strong-coupling regime with large detuning was dis- 
cussed in Ref. [1(| by using a perturbation method . 

Recently, JC model in the intermediate coupling 
regime has been realized in circuit quantum electrody- 
namics (QED) where the superconducting qubits play the 
role of artificial atoms |lll4l5j | . By enhancing the induc- 
tive coupling, the atom-cavity coupling strength reaches 
a considerable fraction (maximum is around 12% ) of the 



cavity transition frequency |13|. Theoretically, some ap- 
proaches are proposed toward the accurate solution to 
the JC model without the RWA flll - [2~Lj ]. Among them, 
the effective approaches within extended bosonic coher- 
ent states have been developed 0, [2l[ . 

In this Letter, we study the JC model without the 
RWA beyond the weak coupling regime by the numeri- 
cally exact solutions [20]. We consider the resonant case 
which is important in the interaction of light and the 
matter. The effect of the CRTs on the dynamic evolu- 
tion of the atomic population inversion and the vacuum 
Rabi splitting is investigated. 

The Hamiltonian of a two-level atom with transition 
frequency uj eg interacting with a single-mode quantized 
cavity of frequency uj is 



H = ~^~ <J z + 0Ja [ a + g(a l + a)a x , 



(1) 



where g is coupling strength, a x and <j z are Pauli spin- 
1/2 operators, oJ and a are the creation and annihilation 
operators for the quantized field. 

By using the extended coherent state technique (20l[2l| . 
we may write a transformed Hamiltonian with a rotation 
around an y axis by an angle ^ 



Hr 



-a x + wa' a + g{a) + a)a z 



(2) 



Then we introduce two new boson operators A = a + g/u 
and B = a — g/u) to eliminate the linear term. The Fock 
states with occupation number m in A and B are 



|m) A(fl) =(at±ff/w)*V 



'|0>, (3) 



The wavefunction of the transformed Hamiltonian can be 
expanded as 



\<f) 



m=0 



|to)a|1) + c m ,2|m)s|2), 



(4) 
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where |1) = ^(|e) + \g)) and |2) = ^(-|e) + with 

|e) and \g) the upper and low atomic state, N tr is the 
truncated number. The wavefunction in original photon 
Fock space can also be expressed as 

\<p) = —, = ' S ^2,{n\{c m ,i\m) a - Cm,2\m) s)\n, e) 

^ n.m 

+ {n\(cm,x\ m ) A + c m M m ) B)\n, g) (5) 

Associated with JC Hamiltonian with CRTS is a con- 
served parity Y[ = e~ l7vN with the excitation number 
N = a^a + \a z + ^. Yl possesses two eigenvalues 
±1, depending on whether the excitation quanta num- 
ber is even or odd. The coefficients therefore satisfy 
c 7rit 2 — ±(— l) m c mj i and the wavefunction with even or 
odd parity can be written as 

- J2 c ™,iHa\1) ± (-IJViHbI^ (6) 

m— 

The coefficients, c m i and the corresponding energies, 
are determined by 

N tr N tr 

E ± Cm,i\m) A = T^f £(-irc ro ,i|m> B | fl > 

m— m— 

iVtr 

+ u^(iU- S 2 /(,; 2 )c roi i|m)i. 

m=0 

(7) 

The solutions can be obtained by digonalizing the set of 
equations with dimension M = N tr + 1. The truncated 
number N tr is set to 40 in the present numerical cal- 
culations, yielding the eigenenergies with relative errors 
less than 10~ 6 . The energies of the first six eigenstates 
and the average photon number in the ground state ver- 
sus g/uj are plotted in Fig. [1] Here and in the following 
calculation, we take the resonant case uj eg = to 

The eigenfunction If^) can be expressed in the states 
|n, e) and n, g) by using the relation (n\m)A = (—l) n D nm 
and (n\rn)B = (—l) m D nm with polynomials D nm = 
S^o" [ "' ml (-l) % !(m ^g- t)! (g/^) m+ "e^ 2 / 2 " 2 , for the 
even parity 

|<^ + ) = ^2 -Cm,iD nm \n,e) + ^2 Cm,iD nm \n, g) 

n—odd,m n—even 

and for the odd parity 
|<fO = c m ,iD nm \n, e) - ^2 c m.\D nm \n,g). 

n—even,m n—odd 

Based on the above new states, Eq.([3]), we study the 
influence of the CRTs on the Rabi splitting, particularly 
the vacuum Rabi splitting. When the atom dressed by 
the cavity mode is subjected to the interaction with the 
vacuum modes, we would have spontaneous emission if 
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FIG. 1: The first six energy levels (black solid line) and the 
mean photons number of the ground states a + a of JC model 
with CRTS (blue squared line) for atom and cavity on reso- 
nance U) = LU eg . 

we pump the dressed atom from its ground to an excited 
state. Under the RWA, when the atom is excited by 
the operator, V = \e)(g\ + \g)(e\, from the ground state 
|<7, >, the emission spectrum has two peaks with equal 
height (the distance of the two peaks, 2g, is the vacuum 
Rabi splitting). When the CRTs are included, the pho- 
tons in the ground state are no longer a vacuum state, 
as shown in Fig. [T] Here, we still use V to excite the 
atom from the ground state of the atom-cavity system, 
and regard the splitting in the emission spectrum as the 
vacuum Rabi splitting. The excited state, V\gs), can be 
expanded in terms of the eigenstates with odd parity 

V\gs) = (\e)(g+\g)(e\M) 

= J2^k)(<Pkl ~c mA D nm \n,g) 

k n—odd,m 

+ ^ c m,iDnm\n,e)]. (8) 

n=even 

From the excited state, the atom will decay to the dressed 
ground state with an emission spectrum. The spectrum 
has Lorenzian peaks due to the spontaneous transitions 
between high energy eigenstate components in the ex- 
cited state to the dressed ground state, and the heights 
of the peaks is proportional to the square of the proba- 
bility of the corresponding eigenstates. The width of the 
peaks is the same determined by the decay rate. The fre- 
quency difference(s) will give us the vacuum Rabi split- 
ting. In Fig. [2] we present the emission spectrum for 
different coupling strength g at resonance u eg =ui. Three 
new features are obviously observed. 

• At small g/oj = 0.1, there are two peaks, but their 
heights are not the same. 

• More than 2 peaks can emerge in the emission spec- 
trum, as shown in Fig. [5] (b) for g/uj — 0.8 with 3 peaks. 
The distance between peaks can be different. The first 
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The atomic population inversion can be expressed as 



g/co =0.1 
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FIG. 2: (Color online) Emission spectrum vr/uj the initial 
excited state V\gs) with to — uj eg for (a) g/uj — 0.1, (b) g/uj = 
0.8. 



splitting is 1.18a; and the second splitting is 0.697a;. 

• At the strong coupling, g/uj = 2, there is no Rabi 
splitting, which is the result of almost equal energy sep- 
aration for all eigenstates. The energy difference of the 
ground state and the first-excited state is 0.0004a; ob- 
tained by Eq.([7]) numerically, and the transition between 
them is smeared out for temperature T > 3 x 10~ 15 o;. 

Without the interaction between atom and the cavity 
mode, we have energy degeneracy for |e, n) and \g, n— 1) 
, and no degeneracy for \g, 0), which results in three peak 
structure for and the vacuum Rabi splitting (two 

peak structure) for n — 0. The degeneracy is broken 
by the interaction. For small g/uj , the non-degeneracy is 
usually kept for all energies under and without the RWA, 
and the eigenenergies without RWA are close to that un- 
der the RWA. For g/u > 0.1, the difference between RWA 
and without RWA merges. For very large g/tu, the degen- 
eracy reappears with different fashion: almost the same 
energy for \<pf) and \ip~). The vacuum Rabi splitting 
disappears. At the same time, the multi-peak structure 
for fi^O also disappears. 

Next, we study the effect of the CRTs on the dynamical 
evolution of the system. Under the RWA we have the 
well-known collapse and revival. It has been shown that 
the collapse and revival would become not clear when the 
CRTs are included @. We investigate the effect of the 
CRTs by increasing g/u with the initial state, the atom 
in the ground state \g) and the field in a coherent state 
| a) with the mean photons n — (a|a^a|a). 

The time evolution of the wavefunction is expressed as 



\<p(t)) =e- iHt \<p(0)) =C7t^;/± e -^%±), (9) 



±k 



where the coefficients are evaluated by / 
(ip^\U\(p(0)). By tracing over the field degrees of the 
freedom, the evolution of the reduced density matrix for 
the atom is given by 



A I 



Pa(t) = l^fkfl 



fc,i 



« ) i Tr pft (C7t|^± )( ^±| t7) . 



P(t) = \Tr((\e)(e\p a (t))\ - \lr((\g)(g\p a (t)). (10) 

with the probability Pj — \(j\p a (t)\j)\ for the atom in 
the upper j = 1 (down j — 2) state at time t. For 
convenience, we also write the population inversion in 
the case of the RWA 



P(t) oc cos(2gWn + l). 



(11) 



Note that the collapse and revival evolution is indepen- 
dent of g in the dimensionless time scale of r = 2gt. 

By using the accurate eigenvalues {E k } and eigen- 
states {<p k } and initial mean photon h — 10, we cal- 
culate the population inversion evolution P(t) without 
the RWA, see the black curves in Fig. [3] for different 
g/oj. For comparison the RWA results are also plotted, 
see the green curves. In the weak coupling regime, say 
g/ui = 0.02, the collapse and revivals can still be seen, 
with some oscillation especially in the collapse period. 
Due to the CRTs, the dependence of 2g^/n in Eq. (11) 
is broken. Therefore, as the increase of g/u> from 0.01 
to 1.0, the collapse and revivals will gradually disappear. 
In Fig. EJb) for g/u — 0.2, rapid oscillation arises in the 
collapse period and the envelope in the revival period 
is smeared, which is consistent with perturbation results 
[9j. As g/u> increases, the effect of the CRTs becomes 
dominated. At g/uj = 0.5, the collapse and revival com- 
pletely disappear due to the CRTs, see Fig. Etc). The 
situation becomes quite different if g/u> further increases. 
For large g/uj = 2, as shown in Fig. 3(d) that the regular 
collapses and revivals reappear, which will be explained 
below. 

We now present the perturbation theory in the strong 
coupling regime. Note that, without the first term 
Hi = —^Y~a x , the Hamiltonian ( [2]) can be diagonal- 
ized exactly in terms of the new operators A and B. The 
kth. eigenstate and eigenvalue read 

l fc± >-^(±(-W B >)'^^ fc ^ 2 /- (12) 

Regarding Hi as a perturbation term, we then obtain the 
second-order perturbative results for energy levels as 



E±=uk-gV^L kk + \^YY. lLlk{2 



V 2 ' fgw(l-k) 



, (13) 



with polynomials L\ k dressed by overlap of two coherent 
states as L kl = (k A \l B )(-l) 1 + (-l) k (k B \l A ). 

By using the above unperturbated eigenstate and 
eigenvalues, the probability of the state with the atom 
in state \i) (i — e,<?) and n photons in the cavity mode 
is given by 



P n ,i{2) = jY.\( anT V™ e ~ M ™~ m) *)l 



(14) 
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FIG. 3: Population inversion evolution 7. (solid curves) for 
the atom-cavity system without RWA for different coupling 
strength g/u = 0.02 (a), 0.2 (b), 0.5 (c) and 2.0 (d). For 
comparison, the corresponding RWA results (green curves), 
and the approximated result from Eq. (|15[) (red curves) are 
also plotted. The initial mean photons n = 10 and cj eg — u. 



It is just the summation of cos[w(n — m)t] that results 
in the recovery of the collapses and revivals, similar to 
the case of the RWA described by Eq.( [IT]). By using 
Eq.f [T5|). the strong coupling perturbative results with 
the red curve is also exhibited in Fig. [3jd) , which almost 
superpose with the black curves. This reappearance of 
the collapses and revivals is originated from the nearly 
equal level spacing for the dressed eigenstates in the 
strong coupling regime, as also demonstrated in Fig. [T] 



In summary, within the bosonic coherent-state ap- 
proach, two interesting features of the JC system well 
studied previously are examined in this paper. Even in 
the present experimentally accessible regime, say, g/uj = 
0.1, the height of the two peaks of the vacuum Rabi split- 
ting are different, the collapses and revivals in the evo- 
lution of the population inversion completely disappear. 
In the deep strong coupling regime, e.g. g=2, the vac- 
uum Rabi splitting vanish and the collapses and revivals 
in the evolution of the atomic population inversion reap- 
pear. Further experiments are then motivated to confirm 
these theoretical predictions. 



with coefficients a n = (nA\a) 2 + (ns|a) 2 ,6„^„ = 
(nA\o>}(riB\a}((rnB\n) a + (mA\n)B)- Therefore, the 
atomic population inversion for strong coupling strength 
is approximately obtained as 

p{t) ~ y _^ a nKm CQs[uj{n _ m)t] (15) 
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